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ABSTRACT 

The recent discovery of terrestrial-type organic species such as methyl formate 
and dimethyl ether in the cold interstellar gas has proved that the formation of 
organic matter in the Galaxy begins at a much earlier stage of star formation 
than was thought before. This discovery represents a challenge for astrochem- 
ical modelers. The abundances of these molecules cannot be explained by the 
previously developed "warm-up" scenario, in which organic molecules are formed 
via diffusive chemistry on surfaces of interstellar grains starting at 30 K, and 
then released to the gas at higher temperatures during later stages of star for- 
mation. In this article, we investigate an alternative scenario in which complex 
organic species are formed via a sequence of gas-phase reactions between precur- 
sor species formed on grain surfaces and then ejected into the gas via efficient 
reactive desorption, a process in which non-thermal desorption occurs as a result 
of conversion of the exothermicity of chemical reactions into the ejection of prod- 
ucts from the surface. The proposed scenario leads to reasonable if somewhat 
mixed results at temperatures as low as 10 K and may be considered as a step 
towards the explanation of abundances of terrestrial-like organic species observed 
during the earliest stages of star formation. 



Subject headings: astrochemistry - ISM - molecular processes 
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Introduction 



Although it has been understood for some time that a combination of ion-molecule 
and neutral-neutral reactions can explain much of the exotic "carbon-chain" chemistry that 
occurs in the gas phase of cold interstellar cores, both starless and prestellar, the chem- 
istry of hot cores, in which gaseous terrestrial (partially saturated) organic species of six 
or more atoms, known as "complex organic molecul es", or COMs for short, are formed 
at relatively high abundance, is not well understood ( jHerbst fc van Dishoeckll2009l ). Well- 
known examples of COMs include methanol (CH3OH), dimethyl ether (CH3OCH3), ethyl 
cyanide (C2H5CN), and methyl formate (HCOOCH3). Early theories of the formation of 
these molecules were based on a warm gas-phase chemistry during the hot co re stage, pos 



sibly preceded by a cold ice-mantle chemistry du r ing the prior cold c ore stage (|Brown et al. 



1988 



Charnlev et al.lll992l ; lHasegawa et al.lll992l ; ICaselli et al.lll993 : Horn et al.ll2004f ). The 



cold ice chemistry is thought to be dominated by reactions involving weakly-bound atoms 
such as hydrogen, which can diffuse at low temperatures. More recently, the chemistry oc- 
curring during the warm-up leading to the formation of a hot core has been emphasized. 
In this approach, the chemical synthesis of hot core molecules begins on the ice mantles of 



grain s in cold cores, where species as complex as methanol are synthesized (IGarrod fc Herbst 



20061 ). As the cold core collapses isothermally into a pre-stellar core and then begins to heat 
up as a protostellar core is formed, the collapsing gas and dust also increase in temperature. 
Starting at temperatures of approximately 30 K, the heavy species formed on the 10 K ice 
begin to diffuse and collide with one another, although they are mainly unreactive. Photons 
or energetic particles strike the dust particles and dissociate some of these heavy species into 
radicals, which are likely to be quite reactive. For example, methanol can be dissociated 
to produce the methoxy radical (CH 3 0) while formaldehyde can form the formyl radical 
(HCO). These two radicals can possibly undergo what is known as a recombination reaction 
to form methyl formate (HCOOCH3). 

Many other COMs are likely formed via recom b inatio n of radicals on ice s duri ng the 
warm-up period, as discussed by IGarrod fc Herbstl ( 120061 ) and IGarrod et al.l (120081 ) . who 
showed that large abundances of COMs in the gas can be achieved by the stage at which the 
temperature of the collapsing material has risen to 100 - 300 K, where thermal evaporation 
takes place. At temperatures above 30 K but below 100 K, smaller abundances of complex 
species formed on grains can also be produced via non-thermal desorption as long as the 
radical-radical recombination reactions can occur. Indeed, observations in various sources 
have shown that small abundances of gaseous COMs can be found in coo ler sources than ho t 
cores such as the galactic center or along the line of sight to protostars (jOberg et al.ll2010l ). 
But, below about 25 K, it would appear that the diffusion of radicals on ice surfa ces can 
no longer occur sufficiently rapidly to react to form COMs ( IGarrod fc Herbstl 120061 ) . Thus, 
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based solely on the warm-up radical mechanism, one would not expect to find any molecules 
of this type in the gas of cold cores such as TMC-1. 

Two recent observations of a number of terrestrial-type complex organic molecules 
i n the cold prestellar cores L 1689b and Bl- b have excited the astrochemical community 



(IBacmann et al.ll2012t ICernicharo et al.ll2012l ) although some more limited work had reached 
similar conclusions for a variety of sources including TMC-1 regarding methanol and ac- 
etaldehyde (CH 3 CHO). New observations of other cold cores (Bacmann, private communica- 
tion(2013)) also show the presence of COMs. The COMs seen in L1689b and/or Bl-b include 
methanol (CH 3 OH), acetaldehyde (CH3CHO), methyl mercaptan (CH 3 SH), dimethyl ether 
(CH3OCH3), and methyl formate (HCOOCH3), while smaller species include the methoxy 
radical (CH 3 0), detected in Bl-b for the first time, the formyl radical (HCO), formic acid 
(HCOOH), ketene (CH 2 CO), and formaldehyde (H 2 CO). Of the complex species, only the 
production of methanol is well understood; it is formed via non-thermal desorption following 
the succ essive hydrogenation by ato mic hydrogen on ice mantles of CO that is accreted from 



the gas (IWatanabe fc Kouchil 120021 ). In this paper, we propose and investigate a scenario for 
the formation of the complex terrestrial-type organic molecules methyl formate, dimethyl 
ether, and acetaldehyde in cold cores. We also report our results for the gas-phase abun- 
dances of ketene, formaldehyde, methoxy, and methanol. The scenario involves an enhanced 
rate of non-thermal desorption of selected precursor molecules (e.g. methoxy, formaldehyde, 
methanol) from cold ice mantles by a process known as reactive desorption, in which the 
exothermicity of surface chemical reactions is at least partially channeled into kinetic energy 
needed to break the bond of the product with the ice. The enhanced desorption rate is 
followed by gas-phase reactions leading to the COMs, among which are radiative association 
processes with rate coefficients that are still poorly understood. Extension of this approach 
to other COMs is certainly feasible. 

The remainder of this paper is organized as follows. In Section 2 we describe our chemical 
model, while in Section 3 modeling results are presented and compared with observational 
data. Section 4 contains a general discussion, while Section 5 is a summary of our work. 



2. Chemical model 

For this study, we utilized a gas-grain chemical mod el with a network of gas-ph ase and 



grain-surface chemical reactions taken from that used by lVasyunin fc Herbstl (120131 ) . which 



in turn is based on the well-established OSU astrochemical databasaH The diffusive surface 



1 http://www. physics. ohio-state.edu/^ eric/research, html 



-4- 



chemistry is treated by the rate equation technique (IHasegawa et al.lll992l ). in which no 
distinction is made between the reactive surface of an ice mantle and the inert icy bulk 
of interstellar grains, which is a no ticeable simplification in c omparison with the structure 

of real multilayer interstellar ices ( Vasyunin fc Herbstl 2013 ). Since the methoxy radical 

(CH 3 0) was detected in the cold core Bl-b ( jCernicharo et al.ll2012l ). we introduced four gas- 
phase reactions for this species, which is produced during the hydrogenation of CO on grain 
mantles followed by desorption. In our network, the methoxy radical and its isomer, the 
hydroxymethyl radical (CH 2 OH), are not distinguished. The first two reactions are simple 
destruction reactions involving the abundant atoms H and O: 



CH3O + H ->■ CH 3 + OH, 



and 



CH3O + O -»■ H 2 CO + OH. 



(2) 



Both of these reactions have been studied to some extent in the laboratory for methoxy and 
hydroxymethyl radicals, as can be seen in the NIST Chemical Kinetics DatabasdE For the 
reaction with atomic hydrogen, only the process involving CH 2 OH has a listed and reviewed 
rate coefficient, which possesses a temperature independent value of 1.6xlCT 10 cm 3 s^ 1 in 
the range 300-2000 K. We assume that the rate coefficient remains constant at temperatures 
down to 10 K in our model. For the reaction with atomic oxygen, measured and reviewed 
values both exist at room temperature for CH 2 OH only; we assume that the average result, 
1.0xl0~ 10 cm 3 s^ 1 , pertains down to 10 K in our model. 

The third and fourth reactions occur between the methoxy and methyl radicals: 

CH3O + CH 3 -> H 2 CO + CH 4 , (3) 



and 



CH3O + CH 3 ->■ CH3OCH3 + hv. (4) 



The rate coefficient of reaction (j3J), also taken from the NIST Chemical Kinetics Database, 
has been measured to be 4.0 xlO -11 cm 3 s _1 by several investigators at room temperature 
and higher. We assume that this value pertains down to 10 K. Reaction (J1J) is a radiative 
association, which has not been studied in the laboratory to the best of our knowledge. The 
system has been studied at high density, where it essentially reaches the collisional limit at 
room temperature. The radiative mechanism, which is the low density process relevant to 
the interstellar medium, is likely t o be much sl ower at room temperature, but to have a 
dependence on temperature of T~ 3 (jHerbstlll980l ) so that by 10 K, reaction (0J is also quite 



2 http://kinetics. nist.gov/kinetics/Search.jsp 
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rapid. Based on previous analyses, we estimate the rate coefficient for this reaction to be 
10~ 15 (T/300 K)~ 3 cm 3 s _1 with a large uncertainty. For convenience, information on new 
and other important reactions is summarized in Table [TJ 

The elemental abundan ces available for gas - grain chemistry in our simulations corre- 
spond to the EAl set from IWakelam fc Herbstl (120081 ). These are so-called "low metal" 
abundances, which are lower by a factor of ~100 for heavy elements in comparison to solar 
elemental abundances in order to take into account elemental depletion on grains in cold 
interstellar cores. The ionization rate ( (s -1 ) due to the main ionizing agent — cosmic rays 
- has a standard value of 1.3xl0~ 17 s" 1 for molecular hydrogen. The initial abundances are 
all atomic except for hydrogen, which starts in the form of H2. The gas-to-dust mass ratio 
is 0.01. 

To treat the diffusive surface chemistry, we assume that the ice mantle surrounding 
interstellar grains can be represented as a two-dimensional lattice with a periodic potential. 
Wells of the potential are the binding sites for atomic and molecular adsorbates. Due to ther- 
mal hopping, species can move across the grain surface and react with each other. In the 
case of atomic and molecular hydrogen, quantum tunneling through the potential barriers is 
also considered as a source of surface mobility. The width of a potential barrier between two 
binding sites is assumed to be 1 A. The height of a barrier against diffusion for a species 
is assumed to be 1/2 of the desorption energy Ed. In our simulations, we consider spherical 
grains of a single size of 10 -5 cm with ~10 6 binding sites. Note that our network is inherited 
from previous studies in which complex organic molecules were formed at higher tempera- 
tures via diffusive surface chemistry. It contains surface radical-radical chemical reactions 
that lead to the formation of COMs. Since at 10 K, the mobility of molecular species on sur- 
faces is negligible, these reactions do not contribute to the formation of COMs in this study. 
The only efficient diffusive reactants on grain surfaces at 10 K are atoms such as hydrogen, 
oxygen and nitrogen, especially atomic hydrogen, which can be an efficient tunneler. For 
example, hydrogen atoms that accrete on to dust particle surfaces are able to hydrogenat e 
CO ice all the way to methanol ice at 10 K (jCharnley et al.lll997l ; IWatanabe fc Kouchill2002l ): 



CO -» HCO HoCO H,CO CH,OH. 



(5) 



The chemistry in the gas phase is connected to grain-surface chemistry via accretion and 
desorption. The accretion probability for all neutral species is set to unity except for helium, 
for which it is taken to be zero because it is chemically inert in surface chemistry and has 
a very low desorption energy, so that it would thermally desorb rapidly anyway. Molecular 
ions are not allowed to stick on grains. There are f our types of desorp tion included in ou r 
model: thermal evaporation ( IHasegawa et al.lll992[ ). photo desorption jFavolle et al.l 120111) . 



desorption due to grain heating by cosmic ray bombardment ( IHasegawa fc Herbst 



1993|), and 
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reactive desorption (IGarrod et al.l 120071 ). In cold dark cloud cores, which are well shielded 
from UV radiation and have very low temperatures of ~10 K, thermal evaporation is nil 
except for He, H, and H 2 while photodesorption can only occur slowly via photons produced 
by energetic electrons. Desorption due to cosmic rays has limited efficiency and has only bee n 
proven to be important on large timescales for species such as CO (IHasegawa fc Herbstlll993l ). 
Constraints on the efficiency of reactive desorption, in which the products of an exothermic 
surface reaction are immediately ejected into the gas, are poor. To the best of our knowledge, 
there are no experimental constraints on its efficiency except for the case of H q dKatz et al. 
1999), and only a simple statistical theory is available for larger products (IGarrod et al. 



20071 ). In this study, we consider three different efficiencies for reactive desorption: 0%, 1%, 
and 10% per exothermic reaction on a grain surface. The 1% efficie ncy is sufficient to e xplain 
observed fractional abundances of methanol (~ 10 -9 ) in co ld cores (IGarrod et al.ll2007ri . The 
10% value is equal to the highest efficiency considered by IGarrod et al.l ( 120071 ). Our models 
with these three different efficiencies for reactive desorption are labeled M0, Ml, and M10. 

All chemical models are run for 10 6 yr, which is slightly longer than so-called "early 
time" , at which the ion-molecule chemistry best reproduces the carbon-chain species observed 
in cold cores. We utilize physical conditions typical for sources such as L 1689b and Bl-b, 
where COMs were found at low temperature: T=10 K, a proton density nn=10 5 cm -3 , and 
a visual extinction Ay=l0. In our analysis, we focus on the time span 10 5 — 10 6 yr. 



3. Results 



Because the calculated ice mantle composition in all three models is quite similar, we 
present the composition for only one model - M10 - in Figure [TJ The similarity of the ice 
composition in all considered models arises because in this work reactive desorption affects 
the rates of all surface reactions by a relatively small factor of at most 10%. As such, it 
does not change considerably the balance between different surface chemical processes, and 
the resulting ice composition remains almost unchanged. As expected, our chemical model 
produce s a typical ice composition which is similar to tha t calculated in a number of previous 
studies ( Garrod fc Pauly 2011 : Vasyunin fc Herbstl 2013), and wh ich is also generally similar 
to observed ice abundances towards protostars (jOberg et al.ll201ll ). The major ice component 
is water, which has a nearly constant fractional abundance of ~10 -4 after 10 5 yr. The 
abundances of the carbon-bearing major ice components change with time more significantly. 
Carbon monoxide accretes from the gas phase and at 10 5 yr possesses an abundance on grain 
surfaces similar to that of water. By 10 6 yr, its abundance has decreased by approximately 
a factor of 3 as it is gradually converted to methanol. Correspondingly, the abundances of 



-7- 



formaldehyde, methanol and methane rise with time by an order of magnitude or more to 
reach values of roughly ten to thirty percent of the abundance of water ice. Another major 
ice constituent observed in ices towards star-forming regions, carbon dioxide, has a very small 
fractional abundance of ~ 4 x in our model, despite a much larger observed abundance. 
This discrepancy arises because carbon dioxide is probably formed durin g other stages of 
star formation when dust grain s have higher temperatures of 
201ll ; IVasvunin fe Herbstll2013h . 



20 K (e.g., [Garrod fc Paulv 



In Figure [21 abundances of selected gas-phase species are presented for models MO, Ml, 
and M10, with formaldehyde, methanol, and methoxy in the left panels and methyl formate, 
ketene (CH 2 CO), dimethyl ether, and acetaldehyde in the right panels. The upper row of 
panels represents results for MO (no reactive desorption), the middle row results for Ml 
(1% reactive desorption efficiency), and the lower row results for model M10 (10% reactive 
desorption effici ency). The re s ults f or Ml are similar to those of the best-fit model in the 
earlier paper of iGarrod et al.l (120071 ). Fractional abundances of all these species observed 
in L1689b and Bl-b are shown in the two rightmost columns of Table [2j We derived the 
observed fracti onal abundances in the se so urces by dividing the o bserved column densities 
as presented in iBacmann et al.l (120121 ) and ICernicharo et al.l (120121 ) by the column densities 
of molecular hydrogen in L1689b and Bl-b. The observed fractional abundances in L1689b 
range from approximately 10~ 10 to 10 -9 whereas those in Bl-b range from approximately 
10~ n to more than 10 -9 . Given the uncertainties both in observed abundances and in our 
astrochemical models, we consider a model to explain successfully the observed abundance 
of a species if the modele d abundance differs from the observed abundance by not more than 
one order of magnitude (jVasyunin et al.l 12004 . 120081 ; IWakelam et al.ll2010l ). 



3.1. MO: a model without reactive desorption 

Formaldehyde is mainly produced by the gas-phase reaction O + CH 3 — > H 2 CO + H, 
while ketene (CH 2 CO) has two paths of formation. First, it is produced in the dissociative 
recombination of C 2 H30 + , which in turn is created in the radiative association reaction 
between the second most abundant gas-phase molecule CO and the simple ion CH^ : 

CO + CH+ ->■ C 2 H 3 0+ + hv. (6) 

Secondly, ketene is produced in the neutral-neutral reaction 

O + C 2 H 3 CH 2 CO + H. (7) 

Other observed organic molecules in the model M0 have very low abundances, because they 
do not have gas-phase routes of formation efficient enough at 10 K to compete with accretion. 
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3.2. Ml and M10: models with reactive desorption 

Our results change drastically when reactive desorption is enabled. For model Ml, as can 
be seen in the left middle row of Figure [21 formaldehyde and methanol reach peak abundances 
in excess of 10 -9 , while the methoxy radical reaches a peak abundance of 10~ 10 . As discussed 
earlier, methanol is not formed in the gas, but is formed by surface hydrogenation of CO 



followed by desorption (IGarrod et al.l 120071 ). This is not only true for CH3OH but also for 
intermediate products of the reaction chain in eq. (jSJ), and other products of surface reactions. 
As such, reactive desorption becomes a major source of methanol and the methoxy radical 
in the gas phase. The abundance of formaldehyde is affected to a lesser extent because 
it is already produced efficiently by gas-phase processes. The increase in the gas-phase 
abundances of methoxy, methanol and formaldehyde after 3x 10 5 yr is caused by the increased 
abundances of these species in the ice mantle, as can be seen in Figure HJ 

For the molecules in the right middle panel, model Ml also shows enhanced abundances 
over MO. While the final abundance of ketene (CH 2 CO) increases only by a factor of three, 
the abundance of acetaldehyde (CH3CHO) is increased by four orders of magnitude, and the 
abundance of dimethyl ether (CH3OCH3) is increased by more than ten orders of magnitude. 
The fractional abundances of both species with respect to hydrogen now exceed 10 -12 towards 
10 6 yr. The additional acetaldehyde is produced via the gas-phase reaction O + C 2 H 5 — > 
CH3CHO + H. The C 2 H 5 product is ejected into the gas-phase via reactive desorption as an 
intermediate product of the chain of surface reactions that hydrogenate C 2 into C 2 He: 

C 2 — > C 2 H —> C 2 H 2 — > C 2 H3 — > C 2 fi4 — y C 2 H5 — > C 2 Hg. (8) 

Gas-phase dimethyl ether, on the other hand, is mainly produced via the gas-phase radiative 
association reaction shown in eq. fl4]). Although inclusion of the "standard" reactive desorp- 
tion with an efficiency of 1% in the model allows us to produce complex organic molecules 
in much higher abundances than with no reactive desorption at all, some of the predicted 
abundances are still quite low. For example, methyl formate (HCOOCH3) has an abundance 
in model Ml below 10~ 14 in the entire time span 10 5 -10 6 yr, which is 3-4 orders of magnitude 
below the observed values. 

The results produced by the model M10, with an efficiency for reactive desorption of 
10%, are presented in the bottom row of Figure [2j In this model, the abundances of all 
depicted species are enhanced in comparison with model Ml and especially with model 
MO. In model M10, the contribution to gas-phase H 2 CO from the chain of surface reactions 
in reaction §5§ via reactive desorption is an order of magnitude larger than the gas-phase 
formation routes, which leads to an increase in the peak fractional abundance of H 2 CO to 
6xl0 -8 . The abundances of methoxy and methanol are increased by an order of magnitude 
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in comparison with model Ml, and reach 8.7x10 10 and 3.4x10 8 , respectively, at their 
maximum. 

As can be seen on the right panel in the bottom row of Figure [2J the abundances of 
ketene and complex organic species are increased in model M10 in comparison with model 
Ml by different factors. The smallest enhancement in abundance, a factor of 7 at 10 6 yr, 
is exhibited by ketene. The order of magnitude increase in the abundance of gas-phase 
methanol throughout the time range shown leads to an increase in the abundance of CH^ 
via the reaction Hj]" + CH 3 OH — > ChLt + H 2 + H 2 . The methyl ion can then react with 
CO in reaction (|6]) to produce C2H 3 + , and subsequently CH 2 CO. An order-of-magnitude 
increase in the abundance of acetaldehyde is caused by the increased abundance of C2H5, 
which is ejected into the gas via efficient reactive desorption from the surface reaction chain 
in reaction (jSj). 

The abundance of dimethyl ether in model M10 is increased by two orders of magnitude 
in comparison with model Ml because the abundances of both gaseous reactants (CH3 and 
CH3O) in reaction (j3J) are increased in model M10 by an order of magnitude. The increase in 
gaseous CH 3 derives from the increase in efficiency of reactive desorption in model M10. In 
turn, CH3 in model M10 is mainly produced in the gas-phase via reaction (TjQ). Another route 
for the formation of CH3OCH3 is via the dissociative recombination of protonated dimethyl 
ether, which is produced in the reaction of methanol with protonated methanol: 

CH3OH+ + CH3OH -> CH3OHCH3 + H 2 (9) 



The strongest abundance enhancement in switching from model Ml to model M10 occurs 
for methyl formate, which increases from 4xl0~ 15 to 1.8xl0 -12 at 10 6 yr. The most efficient 
gas-phase formation route for HCOOCH3 in our model is via dissociative recombination 
of H^COHOCH^ . The principal form ation route of this species occurs via the radiative 
association reaction (IHorn et aJj|2004j ) 



H 2 COH + + H 2 CO ->■ H 2 COHOCH+ + hi/. 



(10) 



The product is not normal protonated methyl formate, which has the structure 140(011)0011^ 
Unfortunately there is a tra nsition state barr ier between the prod uct of react i on (fill and 
protonated methyl formate (IHorn et al.l I2004J ) . We assume, as did IHorn et al.l (120041 ). that 
dissociative recombination of the product of reaction (1101) can also form methyl formate + 
H despite the change in structure required. Reaction (|TU|) is similar to reaction flU]) in the 
sense that its rate depends quadratically on the abundance of gas-phase formaldehyde. Con- 
sequently, the gas-phase abundance of HCOOCH3 in the model M10 becomes comparable 
although still somewhat lower than the observed values at times exceeding 3 x 10 5 yr. If the 
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dissociative recombination reaction does not lead to methyl formate, another process that 
can p roduce protonated methyl formate is the radiative association reaction (IHorn et al. 
20041 ) 

CH+ + HCOOH HC(OH)OCH+ + hi/, (11) 

followed by dissociative recombination. This reaction, however, is not as efficient in produc- 
ing methyl formate in our current models. 

The general be havior of models with reactive desorption is similar to that described in 
Garrod et al.l ( 120071 ). These authors found that at <2-10 6 yr, the abundances of the majority 
of gas-phase species are not affected significantly by the inclusion of reactive desorption. 
The exceptions are several hydrogenated species such as methanol and formaldehyde. In 
this study, we observe a similar picture extended by the enhan ced abundanc e s of C OMs, 
descendants of H2CO and CH3OH, which were not considered in iGarrod et al.l (120071 ). 



3.3. Comparison with observations 



It is interesting to investigate how well our models reproduce the observed abundances 
of complex organic species and their precursors. To quantify the level of agreement, we 
calculated individual confidence parameters Kj for each i-t h spec ies versus time in models 
Ml and M10 according to the prescription in IGarrod et al.l (120071 ) and then took the mean 
of the individual confidence parameters as the criterion of overall agreement. The formula 
for each individual confidence parameter is 



erfc 



niogjXj-log jX, 
V V2a 



obs\ 



(12) 



where erfc is the complementary error function, X, and X° bs are the modeled and observed 
abundances of species i at time t, and the standard deviation a is assumed to equal unity, 
which means that we estimate the uncertainty in observed abundances to be one order 
of magnitude. With this definition, Ki is equal to unity when the observed and model 
abundances are equal, and quickly goes towards a limit of zero with increasing discrepancy 
between these two values. For example, when the difference is one order of magnitude, Ki = 
0.32, and when the difference is two orders of magnitude, Ki = 0.046. Model M0 is excluded 
from this comparison as its results are obviously too far from observational values. 

The molecules shown in Tableware used in the analysis. In addition to these molecules, 
several other simple species were o bserved towards Bl-b and L 1689b such as CO, HCO + , 
N?H+ NH3, CS and HN CQ/HCNO Ifflrano et al.lll999l : lBacmann et al.lliooil : iMarcelino et al. 
20091 ; lOberg et al.ll2010l ). We did not include them in the comparison set for two reasons. 
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First, the abundances of these species come from different sets of observations, and their 
inclusion will make the set heterogeneous. Secondly, the abundances of these species are 
only slightly different in our models with and without reactive desorption. The resulting 
individual and mean confidence parameters for L 1689b and Bl-b are shown in Figures El and 
HI respectively. For L1689b, model Ml reaches its best agreement at 1.3xl0 5 yr and model 
M10 at 5.1 x 10 5 yr (if we exclude for M10 a small cusp at much earlier time that has a similar 
extent of agreement). In the case of Bl-b, the time of best agreement time for model Ml is 
1.0 x 10 6 yr, while for model M10 it is 2.6 xlO 5 yr. At these times, the average agreement 
is better than one order of magnitude. The individual molecular abundances at these best 
times are shown in Table [2j Because there are two more organic species detected in Bl-b 
than in L1689b, we choose Bl-b for detailed discussion. 

At the time of best agreement for Ml, the modeled abundances for the two precur- 
sor species CH2CO and H2CO, and the two complex species CH3OH and CH3CHO, out 
of seven species, differ from the observational values by less than one order of magnitude. 
For one species, CH3OCH3, the difference between modeled and observed abundances just 
slightly exceeds one order of magnitude. At the same time, model Ml somewhat overpro- 
duces methoxy and significantly underproduces methyl formate. Model M10 exhibits its best 
agreement with observations of Bl-b significantly earlier. While it manages to reproduce the 
abundances of all complex organic species including HCOOCH3, the modeled abundance of 
H 2 CO is two orders of magnitude higher than the observed value. The abundance of CH 3 
is also significantly higher than the observed value. 

A possible reason for the overl y high calculated abund ance of methoxy in both the Ml 



and M10 models has been given by lCernicharo et al.l (120121 ). who showed that that on grains 



the isomer CH 2 OH is produced more efficiently than CH 3 and released into the gas phase. 
This isomer has different spectral properties, and was not observed in the gas phase. The 
authors proposed an alternative formation route for methoxy via the gas-phase reaction 

CH3OH + OH -»■ CH 3 + H 2 0. (13) 

Although our chemical model currently does not discriminate between CH 3 and CH 2 OH, 
we estimate that if methoxy is only produced by reaction ([TBI at its newly determined 
rate coefficient, then the abundance of methoxy will be in much closer agreement with 
observations, at least in case of Bl-b. But, one must remember that since methoxy is a 
precursor species to complex organic molecules, lessening its abundance will also lessen the 
abundance of the complex molecules formed from it, especially dimethyl ether. 

We show the organic abundances vs time for an altered model M10, in which the gas- 
phase route to methoxy is included and the reactive desorption route is removed, in Figure [5j 
The results can be compared with the abundances obtained with the original model M10, 
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which are shown in the bottom panel of Figure [21 The time of the best agreement between 
observed values for Bl-b and results of the altered model M10 is the same as for the original 
model, but at later times, the level of agreement remains at almost the same value until 10 6 yr. 
For Bl-b, the new fractional abundance of CH 3 at the time of best agreement is 5xl0~ n , 
which is closer to the observed value, although still too large by an order of magnitude. This 
fact may imply that the rate coefficient of reaction ([131) is somewhat overestimated, which is 
supported by older estimations. With the high value for the rate coefficient of reaction (T13|) . 
the new abundance of CH3OCH3 at the moment of best agreement for Bl-b is 1 x 10~ 12 , 
but it increases steadily, reaching a value of 2xl0 -11 at 10 6 yr. 

Finally, it is interesting to note that in Bl-b the observed abundances of all organic 
species we consider are lower than in L1689b, including formaldehyde. If future observations 
confirm that this correlation is not caused by systematic errors between observational data 
sets, the lower abundances in Bl-b argue in favor of our scenario for the formation of complex 
molecules in cold clouds, because our scenario implies a strong correlation between the 
abundances of H2CO and COMs. 



4. Discussion 

We have considered one approach to understanding how terrestrial-type organic molecules 
can be synthesized in cold starless and prestellar cores. In our approach, reactive desorption 
succeeds in removing precursor molecules efficiently from the ice mantles of cold dust grains 
into the gas, where subsequent reactions produce the organic molecules in reasonable, if 
not perfect, agreement with the small observed abundances detected in the cold prestellar 
cores L1689b and Bl-b. Radiative association reactions in the gas play an important role 
in the synthesis of the complex organic molecules methyl formate and dimethyl ether. Both 
reactive desorption and radiative association reactions are still poorly understood processes 
in general. Moreover, the product for the proposed radiative association reaction that leads 
to the formation of protonated methyl formate most likely yields a different structure than 
desired, which may or may not form methyl formate in its lowest conformer when the ion 
undergoes a dissociative recombination with electrons. 

Although we have discussed our model results mainly for observed COMs, our model can 
also be tested by its predictions for the abundances of other complex species, such as ethyl 
cyanide (C2H5CN), which is a well-known weed in hot cores. Under cold core conditions, our 
model first produces ethyl cyanide in the ice mantle by hydrogenation of the CCCN radical 



(IHasegawa et al.lll992l ; ICaseUi et al.lll993l ). Reactive desorption then results in a considerable 



abundance of this species in the gas phase at selected times. For example, in Model Ml, 
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the fractional abundance of ethyl cyanide peaks at over 10~ 10 at a time near 2 x 10 5 yr, 
while for Model M10, the peak value approaches a fractional abundance of 10~ 8 at the same 
time. At the times of best agreement for models Ml and M10 with Bl-b, the calculated 
fractional abundances of ethyl cyanide are somewhat under 10~ 12 and somewhat over 10~ 10 
respectively, while for L1689b, these numbers are somewhat over 10~ 10 and near 10~ n . 

There are perhaps other complementary approaches to the production of COMs in 
cold cores, where the dust surfaces are typically to o cold for the surface radical-radical 
photochemistry suggested by lGarrod fc Herbstl (120061 ) . One possibility is that we are missing 
a large number of thermal surface/ice reactions that can occur at 10 K, typically by the 
diffusion of atoms somewhat heavier than hydrogen. A figure showing the many species that 
might be synthesized i n this way was produced by S. C harnley (2009, private communication) 
and can be found in iHerbst fc van Dishoeckl (120091 ). Another possibility is that the dust 
particles are not always at 10 K, and that a variety of effects cause a temporary increase in 
temperature high enough to allow surface radicals to diffuse and react with one another. It is 
even possible that photodissociation or particle bombardment processes can cause temporary 
non-thermal motions of species so as to promot e the rates of non-thermal diffusive reactions 
( jOcchiogrosso et al.l 1201 It iKim fc Kaiser! 120121 ) . Finally, tunneling from one lattice site to 
another for species other than atomic hydrogen will promote more rapid surface/ice reactions 
at lower temperatures than if the only diffusive mechanism for species heaver than H is simple 
thermal hopping. We hope to investigate some of these possibilities in the near future. 



5. Summary 

In this study, we have proposed and examined a scenario of formation of complex organic 
molecules recently found in the cold cores L1689b and Bl-b with gas temperatures of 10 K. 
The proposed scenario is based on the assumption that complex organic species in the cold 
regions are formed via gas-phase ion-molecular and neutral-neutral chemistry from simpler 
precursors such as formaldehyde and methanol. These precursors are formed on icy surfaces 
of interstellar grains, and then ejected into the gas via reactive desorption. Because of our 
lack of knowledge of the efficiency of reactive desorption, we considered three models: M0, 
Ml, and M10, in which the efficiency of reactive desorption per exothermic surface reaction 
is 0%, 1%, and 10%, respectively. 

We found that the proposed scenario gives somewhat mixed agreement with observations 
of COMs and their precursors. Both models with non-zero reactive desorption exhibit much 
better agreement with observations than the model without reactive desorption. The best 
agreement for the whole set of species is exhibited by model M10, which possesses the 



- 14 - 



highest efficiency of reactive desorption. On the other hand, this model tends to overestimate 
gas-phase abundances of formaldehyde in comparison to observations by ~ two orders of 
magnitude. Such a discrepancy can be partly but not fully explained by uncertainties in 
our chemical model, particularly in the rate coefficients of chemical reactions. This fact 
may imply that the proposed scenario of formation of complex organic species via gas-phase 
reactions driven by reactive desorption is only partially responsible for the formation of 
these molecules, and should be assisted by other routes such as processes that warm up 
grains temporarily or that permit non-thermal ice reactive processes. 
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Table 1: New and Other Important Gas-phase Chemical Reactions 
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Notes: a(-b) stands for a x 10 b . a, (3 and 7 are the coefficients in the modified Arrhenius expression for 
the reaction rate coefficient: k = a ■ (T/300)^ • exp(— 7/T). URL for NIST Chemical Kinetics Database: 



http:/ /kinetics. nist.gov/kinetics/index.jsp 
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Table 2: Observed and best-fit modeled fractional abundances of organic species detected in 
cold interstellar clouds 
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Fig. 1. — Abundances vs. time of water and major carbon-bearing ice compounds in the 
model M10. The prefix "gr" denotes species on and in ice mantles of interstellar grains. 
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Fig. 2. — Gas-phase abundances vs time for simple and complex organic species in models 
M0 (top row), Ml (middle row), and M10 (bottom row). 
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Fig. 3. — Mean and individual confidence levels for the species in models Ml (upper panels) 
and M10 (lower panels) in the case of L1689b. The higher the mean confidence level, the 
better the agreement between modeled and observed abundances. 
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Fig. 4. — Mean and individual confidence levels for the species in models Ml (upper panels) 
and M10 (lower panels) in the case of Bl-b. The higher the mean confidence level, the better 
the agreement between modeled and observed abundances. 
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Fig. 5. — Abundances vs time of complex organic species and their precursors in the altered 
model M10, with reactive desorption of methoxy radicals switched off and reaction (fl3|) 
enabled. 



